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Abstract 

We study the evolution of cosmic strings taking into account the 
pJIj ■ frictional force due to the surrounding radiation. We consider small 

perturbations on straight strings, oscillation of circular loops and small 
perturbations on circular loops. For straight strings, friction expo- 
j^j ■ nentially suppresses perturbations whose co-moving scale crosses the 

horizon before cosmological time ~ u~ 2 (in Planck units), where 
fi is the string tension. Loops with size much smaller than t* will 
be approximately circular at the time when they start the relativis- 
tic collapse. We investigate the possibility that such loops will form 
black holes. We find that the number of black holes which are formed 
through this process is well bellow present observational limits, so this 
does not give any lower or upper bounds on u. We also consider the 
case of straight strings attached to walls and circular holes that can 
spontaneously nucleate on metastable domain walls. 



1 Introduction 



Cosmic strings are topological defects that may have formed during phase 
transitions in the early universe (see e.g. Jl||). Their properties and obser- 
vational consequences, especially in connection with their possible role in 
the formation of large scale structure in the universe, have been extensively 
studied during the past decade (see e.g. 0). 

Cosmic strings of mass per unit length \i would have formed at cosmo- 
logical time of order t ~ (Gfj,)~Hpi (G is Newton's constant and t pi is the 
Planck time). It is well known that immediately after the phase transition 
the dynamics of strings would be dominated by the force of friction [^|, ffl. 
This force is due to the scattering of thermal particles off the string. Fric- 
tion would dominate the dynamics until a time of order t* ~ (G/j,)~Hpi. In 
most of the investigations about cosmic strings, the effects of friction have 
been neglected. The reason is that if cosmic strings have to play a role in 
galaxy formation, then they have to form near the grand unification scale. 
In that case their mass per unit length is of order Gfi ~ 10 -6 and friction is 
important only for a very short period of time. 

However, if strings have formed at later phase transitions, say, closer to 
the electroweak scale, their dynamics would be dominated by friction through 
most of the thermal history of the universe. It is therefore of some interest to 
study the evolution of cosmic strings with friction in a quantitative way. The 
relativistic equation of motion for strings with friction was given by Vilenkin 
[[§]. The main purpose of this paper is to solve this equation in a few simple 
cases, which should be representative of more complicated situations. 

The plan of the paper is the following. In Section 2 we review the re- 
sults of Ref. [f|], in order to fix the notation. In Section 3 we study linearized 
perturbations on an infinite straight string. In Section 4 we consider the dy- 
namics of oscillating circular loops. The evolution of linearized perturbations 
on circular loops is discussed in Section 5. 

It is known that an exactly circular loop would form a black hole when it 
collapses under its tension || [7| . Since friction tends to erase perturbations 
whose co-moving scale crosses the horizon before t*, loops smaller than £* will 
tend to be approximately circular. In Section 6 we study the possibility that 
such loops form black holes. Note that these black holes can have masses of 
order M ~ /it*, which can be considerably large if G/i is sufficiently small. 
We discuss possible observational consequences of these black holes. 
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Finally, in Section 7, we consider perturbations on strings which are at- 
tached to domain walls, and in particular, to circular holes which can spon- 
taneously nucleate in a metastable wall. Our conclusions are summarized in 
Section 8. 



2 Strings with friction 

In this section we summarize the results of Ref. |3J, where the equation of 
motion for a string with friction was found. A cosmic string can be repre- 
sented as a two dimensional worldsheet in spacetime x^ = x /i (^ a ), where £ a , 
(a = 0, 1) are two arbitrary parameters on the worldsheet and x^ are space- 
time coordinates (// = 0, ...,4). In flat space, the equation of motion for the 
strings is fiOx u = 0, where □ denotes the d'Alembertian on the worldsheet. 
When one includes the effects of curvature of the spacetime and a force of 
friction F u , the equation of motion reads H 



-□^ + C T x> r,a =F>i,T,a) (1) 



Here T u aT are the four dimensional Christoffel symbols (greek indices run 
over spacetime coordinates). The force of friction F v will depend on the 
temperature of the surrounding matter, T, the velocity of the fluid transverse 
to the worldsheet u\_ = u u — x u a x a,a u a , and the type of interaction between 
the particles and the string, which we symbolically represent by a. Vilenkin 
UD found the form of F v for the case when friction is dominated by Aharonov- 
Bohm scattering of charged particles with the pure gauge field outside the 
string || 

F v = f3T\ u ± . (2) 
The numerical coefficient (3 is given by 

/ 3 = 27r- 2 £(3)5> a sin 2 (7rz/ a ) 

a 

where the summation is over all effectively massles degrees of freedom (i.e. 
mass m « T) b a = 1 for bosons, b a = 3/4 for fermions and v a is the phase 
change experienced by a particle as it is transported around the string. 
For the case of a Friedman-Robertson- Walker (FRW) universe, 

ds 2 = a 2 (r)[-dT 2 + dx 2 ], (3) 
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with the four-velocity of the fluid given by u u = (a 1 , 0,0,0), and choosing 
the gauge £° = r, x ■ x' = 0, Eq.((l|) reduces to 



x' \ r „ , TT , m 1.2s 



f-e" 1 — ] +[2a(H + h)](l-x)x = Q, (4) 



•2 



e + [2a(H + h)]x e = 0. (5) 
Here ' = d/dr, ' = d/d^ 1 , H = a/a 2 is the Hubble parameter and 



2/i 



The quantity e = [a;' 2 /(l — x )] 1//2 is related to the energy of the string [|TJ 

E = /ia(r) y d^e. (6) 

[Using the definition of e it can be easily seen that Eq. (|^) is actually not 
independent of (|)]. 

As noted in Ref.Q, the effect of friction is to make the replacement 

H ^ H + h 

in the equations of motion for a free string in an expanding universe. In what 
follows we will consider 

a(r) ~ r a 

(a = 1 in the radiation dominated era and a = 2 in the matter era). Since 
T ~ a -1 , the friction term h(r) dominates at early times, while the expan- 
sion term H(t) will dominate at late times. It is clear that friction will be 
unimportant for r >> r*, where r* is the time at which both terms are equal 

h( n ) = H(n). 

In terms of cosmological time t defined by dt = a{r)dr and using Ein- 
stein's equations in the radiation era, H 2 = (87r 3 G/90)A/'T 4 , one finds [|J 

U = A(Gfi)-%, (7) 
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(for G/i~(At p /t eg ) 1 / 2 ). Here t eq is the time of equal matter and radiation 
densities and A = [90/? 4//3 /327r 3 A/"] 3//2 , with N the effective number of mass- 
less degrees of freedom. The coefficient A can be rather small. Taking 
sin 2 7iu a ~ 1/2 in the expression for (3 one has 

A ~ 4 x 10" 4 AT 1/2 . 



This estimate should actually be taken as an upper bound for A, since some 
of the particle species may not interact with the string. 

If the strings are so light that Gn^(At p /t eq ) 1 ^ 2 then £* > t eq and (0) is 
not valid. Then we have to take into account that T ~ t~ 2 / 3 in the matter 
era and we have ^ 

UnAWiGp)- 1 ^ V (8) 



3 Linearized perturbations on a straight string 



Consider a straight string at rest, with trajectory x = (£, 0, 0) and introduce 
a small displacement in the plane transverse to the string, 

Substituting in (B) and keeping only terms linear in Sx we have 



y + 2a 



2a- 1 



-2a 



y + k 2 y = 0, 



(9) 



where, as mentioned before, a = 1 in the radiation era and a = 2 in the 
matter era. 

The case without friction corresponds to taking r* = 0. Then (|9]) is 
just the Modified Bessel equation. The corresponding solution that is well 
behaved as r — > is given by 

y = y 2»T{v + l){kT)-»J v {kT), 

where v — a — 1/2. At early times, r << k~ l , the solution tends to a 
constant, 

y^y , (r«k- r ) (10) 



5 



whereas at late times, t » k 1 , the solution is oscillatory 



1)!! . . cos(£:r ). (r»k i ) 

a r) 2 



(11) 



Modes with wavelength larger than the horizon size (r << A;^ 1 ) are frozen in, 
and their physical amplitude y p h ys = o(r)y is conformally stretched by the 
expansion. Once the wavelength of a mode falls within the horizon (r >> 
A;" 1 ), the perturbation starts oscillating with constant physical amplitude |]]. 
When friction is included (r* 7^ 0), equation (g) can only be solved in the 



high frequency approximation {k » 



This is actually the interesting 



regime, in which friction plays a role. In the opposite limit (k << r" 1 ), we 
just saw that perturbations are frozen in up to a time r ~ k~ x » r* due 
to the expansion term 2a /t alone, and by the time they start oscillating the 
friction term has become negligible. 

To solve for the case k >> r~ l we introduce the variable 



r a exp 



a 



2a-l 



y- 



(12) 



Substituting in 



with potential 



2a - 1 

one finds that \l/ satisfies the Schrodinger equation 

- ^ + V(t)V = k 2 V, (13) 



( ( 1 — 1 ) -i ( 1 ( — 



Aa-2 



Now Eq.flTB]) can be solved in the WKB approximation. This is done in 
Appendix A. 

With the boundary condition that the perturbations are frozen in with 
the expansion at early times, the solution is 



y{ r ) 



yo 



y(r) w 2y Q 



a(r k ) 



W(t/t,) -~/kT k 



or 



J « T k ) 



cos(kr + 



(14) 

(r » r k ). (15) 



Here 



a 



2a 



2a-l 
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7 = B 



l + 2a 1 



4a -2 V 4a ' 2, 
where B is Euler's Beta function and Tk is the classical turning point of 
the corresponding Schrodinger problem V(rk) = k 2 . The phase <pk is just a 
constant. In radiation dominated universe 

n = (ArV*) 1 / 2 , 

7 w 1.2 

In general 

r k « [afc-V^- 1 ] 172 "- (16) 

It can be checked that the conditions for the validity of the WKB approxi- 
mation are satisfied provided that k >> r^ 1 and a > 1/2. 

Comparing (|I3) with ( |Tl"D we can summarize as follows. Without friction 
perturbations of wave number k are conformally stretched up to a time r ~ 
A; -1 , after which they start oscillating with constant physical amplitude 

y P hys « (2a - l)!!a(A;- 1 )yo. (17) 

With friction included, short wavelength perturbations k » r^ 1 are 
conformally stretched up to a time Tk given by (|i~6|), hence the factor a{jk) 
in fll5|). Since > A; -1 , friction contributes to increase the amplitude of the 
perturbation by a factor (rkk) a , from time r ~ up to time Tk- After 
Tfc, the perturbation starts oscillating and losing energy to friction. This is 
represented by the two exponential factors in fllip . By the time r ~ the 
amplitude of the perturbations has been damped by a factor exp(— jkTk). 
For r >> r* the modes do not lose any more energy to friction, the first 
exponential in fll^D has reached its asymptotic value of unity and the string 
oscillates with constant physical amplitude 

y phys w 2 ■ a{r k )e-^y . (18) 

As expected, friction plays no role after r*. Also, the exponential sup- 
pression in (|18~D becomes less and less dramatic as approaches r* (i.e. 
A; -1 ~ Tfc ~ r*), in agreement with the intuitive expectation that friction 
does not affect wavelengths of co-moving size comparable or larger than r*. 

As mentioned before, the WKB approximation is only valid for a > 1/2. 
The case a < 1/2 is markedly different because friction never 'switches off'. 
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As a result, the amplitude of perturbations is eventually damped to zero 
for all wavelengths. The simplest example is the flat space case, (H = 0, 
h = constant.) In this case the equation for y is that of a damped harmonic 
oscillator and all excitations disappear after a characteristic lifetime which is 
given by r ~ h^ 1 for k >> h, and by r ~ 2hk~ 2 for k «h (taking a = 1). 



4 Circular loops 

A circular loop can be parametrized as 

x(t,0 = R(t) ■ (cos(e/i?o),sin(e//? ),0), (19) 

where Ro is a constant and £ G [0,2tcR ]. The equations of motion (Q) and 
(H) read (see Appendix B) 

R + 2a(H + h)(l- R 2 )R + ^ ~ ^ = 0, (20) 

R 

e + 2a(H + h)R 2 e = 0, (21) 

with 

R 

(22) 



Ro(l -R 2 ) x / 2 

In the case of flat space and no friction, if = /i = 0,eisa constant, so 
Eq.(^) is linear and one finds the well known oscillatory solution 

R(t) = R cos(t/R ). (23) 

Note that, as a consequence of exact circular symmetry, this solution collapses 
to a point. Of course the same will happen when we include friction and 
expansion. In reality, when the string shrinks to a point it would form a 
black hole 0, [llj. We will consider this process in more detail in Section ^|, 
but for now we shall ignore it, and continue the solutions beyond the singular 
points. The reason is that here we are interested in the energy loss of a loop 
due to friction. Although we consider the circular loop for simplicity, similar 
results would apply to nearly circular loops which do not shrink to a point. 

In the generic case H + h ^ we cannot find exact analytic solutions, so 
we have to use numerical solutions and analytic approximations. 
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a-Loops with friction but no expansion 

It is instructive to start by considering the case H = and h = const. ^ 0. 
Consider a loop initially at rest and whose initial radius Rq is sufficiently 
large. It is clear that as long as aR » h^ 1 the motion of the loop will be 
overdamped, with characteristic velocity 

R « (2ahR)' 1 

[i.e. we neglect the R term in (ppf)]. That this is a good approximation can 
be checked by computing R from the previous equation and comparing it 
to the other terms in (|2l|) . On the contrary, for aR << h~ l the damping 
term can be neglected and the string undergoes a relativistic collapse similar 
to the frictionless one. This suggests that the energy of the string at the 
moment of its first collapse will be independent of Rq (for R » h^ 1 ), and 
will be roughly equal to the energy of a loop of physical size ft, -1 , 

E ~ /i/i" 1 . 

This estimate can be made rigorous by using a simple scaling argument. 
Introducing dimensionless variables p = R, u = 2haR and v = 2haR Q e, the 



differential equations (p0|) and (|2T|) can be reduced to 

+ (l-p 2 



dp 
du 



'1-P 2 ) 1/2 , 2 , 



pv 
dy 



-pv, 



du 

where there is no reference to time or to h. The initial conditions R = Ro 
and R = become 

v(u ) = u , p(u ) = 0, 

where Uq = 2haR Q . The energy of the string at the moment of first collapse 
is 

E = 2ir{iaR e = \ • tc^ 1 . (24) 

The coefficient x = v ( u — 0) can De found by numerical integration. Al- 
though v(u — 0) depends on w , the result rapidly saturates to a constant 
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for Mo larger than 1. 

After the first collapse the loop will undergo a series of oscillations, losing 
energy to friction in each one of them. Let us estimate this energy loss. Since 
the behaviour of energy is controlled by Eq. (|2lD we will be interested in the 
quantity < R 2 >, where the brackets denote the temporal average between 
two consecutive collapses of the loop. We have 

< R 2 >= - < RR >=< (1 - R 2 )[l + 2ahRR] >= 

1- < R 2 > -2ah < RR 3 >, 
or < R 2 >= (1 — 2ah < RR 3 >)/2, where in the first step we have inte- 



grated by parts and in the second we have used the equation of motion (|20|) . 
Repeating similar steps for the calculation of < RR 3 > one can generate a 
perturbative expansion in powers of ahR. The first terms are 

< R 2 >= I _ l -a 2 h 2 <R 2 > +0(a 4 h 4 R A ). 

2 8 V 1 



From (|24] ) it is clear that after the first collapse a/ii?~.3, so the second term 



is a very small correction (at most of order 10 2 ) to the first and we have 

< R 2 >^ 1/2. 



From (^T|) the fraction of the energy lost in between two consecutive col- 
lapses is given by AE/E « —2 < R 2 > hT, where T is the time spent in 
the oscillation. Therefore E oc exp(-ht) and the loop will disappear in a 
characteristic lifetime of order ~ h^ 1 after the first collapse. 



b- Loops with friction and expansion 

The case without friction in expanding universe has been previously studied 
in the literature |L|. Loops are conformally stretched by the expansion until 
they cross the horizon. This happens at time t c defined by 

H-\t c ) = a(t c )R = r c . 

Here r c is the physical radius of the loop at horizon crossing. After that, the 
loop oscillates with constant physical amplitude. This behaviour is illustrated 
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in Fig la, which corresponds to the numerical solution of ( |20| ) with h = and 
for a radiation dominated universe (a = 1). The result is plotted in terms of 
dimensionless cosmological time t/t c . The upper line represents E/E c , where 
E is the energy of the loop and E c = 2irfir c . In the same figure, we plot r/r c , 
where r = aR is the physical radius, and R, the velocity of the string with 
respect to the cosmological fluid. 

With friction included, for loops such that r c >> £*, the evolution is not 
much different from the one just described. The reason is that by the time 
these loops cross the horizon and start oscillating, we already have h « H. 

The effect will be important for loops that cross the horizon well before 
i* (i.e. r c « t*). In Fig. lb we plot the time evolution of a loop with 
r c = 10~ 3 t*, in radiation dominated era (same conventions as in Fig. la). 
Not surprisingly, friction delays the time at which the loop first collapses. 
Initially, this increases the energy of the loop, since it is stretched up to 
a radius much larger than t c . Later, as the loop shrinks, it loses energy to 
friction. It is interesting to observe that both effects roughly compensate each 
other, in the sense that at the time of first collapse, tf we have E(tf) ~ E c , 
just like in the frictionless case. After tf the loop keeps losing energy during 
each oscillation, up to a time i*, when friction switches off. From that time 
on, the loop oscillates with constant physical amplitude. 

It is possible to give an approximate analytical description of this evolu- 
tion. At early times the motion of the loop is overdamped. Neglecting the 
R term in (|2~0|) we have 

RR « ~} — . (25) 
2a(H + h) v ; 



Before £*, H can be neglected and we can integrate (|25|) to find 

r 2a+l 

The loop will first collapse at a value of conformal time which is of order 

r^r f = [a(2a + l)rf ~ l Rl} l /( 2a+1 \ (27) 

Actually, Eq. (|26|) is only valid as long as R « 1. After that, the R term 
in (pOD becomes comparable to the others. Since R$ « Tf, the effective 
friction coefficient, a/t, will not appreciably change during the relativistic 
collapse (which occurs at a conformal time close to Tf). From fl2l|), the 
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energy of the loop at the moment of first collapse is Ef =« .57np/h(rf). 
This has to be compared to E c = 27r//r c ~ 2npR a(Ro). Using fl27l) for 77 
we have 

El w (_ 
E c \Ro, 

In the radiation era (a = 1), we have Ef ~ £7 C independently of Rq, in 
good agreement with the numerical results. In the matter dominated era, 
for R << r* we have Ef » E c , so by the time r/ loops have actually more 
energy than they would have had in the absence of friction. 

After tf, the radius of the loop is much smaller than the horizon and the 
quantity 2a(H + h) is slowly varying compared with the period of oscillation. 
From the arguments of the previous subsection we have < R 2 >~ 1/2, and 
from pip the change in e during one oscillation is Ae/e = —a(H + ft,) At. 
Since Ae << e, one can take the continuous limit: 

e a ar 2 ^ 1 



e a T 2a 



and since E oc ae, the energy of the loop is given by 



E w Ef exp 



t 2q-1 

a^r; — dr 

Tf T 



(28) 



In the case without friction (r* = 0) we have E ~ constant, in agreement 
with the numerical results (Fig. la). With friction included, the energy 
drops until time r#, after which loops oscillate with constant energy. 



5 Perturbations on circular loops 

Let R{t) be a solution of (p0|), representing the evolution of a circular loop. 
A perturbed loop can be parametrized as 

p = R(r)+y^T,9), (29) 

z = y\r,Q). 

Here (p, 9, z) are co-moving cylindrical coordinates, in which the metric takes 
the form ds 2 = a 2 {r){—dr 2 + dp 2 + p 2 dd 2 + dz 2 ), y p is a radial perturbation 
and y z is a perturbation transverse to the plane of the loop. 
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It is straightforward to write (01) in cylindrical coordinates and then sub- 
stitute fl29|) to find the linearized equations for the perturbations. This is 
done in Appendix B. Decomposing y p as a sum over modes 

oo 

y p = E [v p l + Mm + y p f cos(l#)] (30) 

L=2 

and similarly for y z , the resulting equations are 

f L + 2a(H + h)(l- 3R 2 )y p L - + ^^V P L = 0, (31) 

y z L + 2a(H + h)(l- R 2 )y z L + -^y z L = 0. (32) 

where we have omitted the index + or — . The sum in (|30D starts at L = 2 
rather than at L = because it is easy to see that to linear order the L = 
and L = 1 modes do not correspond to deformations from the circular shape 
but rather to small translations and rotations of a circular loop ||. 

Since the function R(t) is not known analytically, Eqs.(j3lD and (|32|) have 
to be solved numerically. However, before we do that, the behaviour of the 
perturbations can be guessed from the results of the previous sections. 

Let us first consider the case without friction. Up to the time t c the 
loop is conformally stretched by the expansion and the perturbations will 
approximately behave as perturbations on a straight string. Ignoring oscil- 
latory cosine factors, the amplitude of the perturbations at this time will be 
yL^c) ~ yo,Ld(Ro/L)/a(r), where y Q L is the initial perturbation. Here we 
have used ([□]) making the substitution k — > L/Rq. After t c , the loop comes 
within the cosmological horizon and starts collapsing. A loop on subhorizon 
scales behaves approximately as it would behave in flat space. The theory of 
perturbations on a circular loop collapsing in flat space was solved in || (see 
Appendix A of that reference). It was shown that the amplitude of transverse 
perturbations stays constant during the collapse, whereas the amplitude of 
radial perturbations shrinks by a factor of L as the loop shrinks to r << r c . 
Therefore by that time we have 

y p phys , L - a(r)vZ - y p o, L ^^-, (33) 

y z ph ys,L = a(r)y z L w y* L a(R /L). (34) 
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To check the validity of these approximations one has to numerically solve 
Eqs.(j3lD and (p2|), whith h = 0. We have done that for different values of L, 
in the radiation dominated universe and using the boundary condition that 
perturbations are initially at rest. The result is plotted in Fig. 2 as a function 
of wave number L. The circles denote the ratio y P vhys L divided by the r.h.s. of 
equation (|33|) at the time when the loop first collapses. Similarly, the crosses 
denote the ratio of y z phys L to the r.h.s. of ([34|). Ignoring the "valleys" in Fig. 
2, these ratios are of order 1, which means that Eqs. (|33|) and (|34]) give a 
very good estimate. The valleys in Fig. 2 are due to oscillatory behaviour of 
the perturbations, which we have ignored in our argument. 

The effects of friction can be introduced along similar lines, and they 
will only be important for co- moving wavelengths (Rq/L) < r*. For such 
wavelengths the r.h.s. of equations ([33|) and ([34]) has to be corrected by a 
factor of 

2 

— — -— {kT k ) a exp(-^kr k ). (35) 

This is obtained comparing ( |17D with where now k is given by L/Rq. 
Again, one can check the accuracy of this approximation by solving the 
equations of motion for the perturbations, now with ft ^ 0. The results 
for the case of a radiation dominated universe are plotted in Fig. 3a for 
Rq = , and in Fig. 3b for R = 5t*. The circles and crosses denote the 
same quantities as in Fig. 2. It is seen that the suppression factor (|35|). 
depicted as a solid line, gives the right answer to very good approximation. 

The above considerations apply only to loops with i? ~' r *- For R << r* 
the perturbations do not have time to be damped before the loop starts 
shrinking, so the factor (^) actually overestimates the effect of friction. We 
shall consider this in more detail in the next section. 

6 Black hole formation 

Let r c be the physical radius of a circular loop at horizon crossing. The mass 
of this loop is 27rr c /i and the Schwarzschild radius corresponding to this mass 
is 

r s = 4irG/ir c . 

As the loop shrinks under its tension, its rest mass is converted into kinetic 
energy so that the total energy of the loop remains constant (neglecting 
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friction and gravitational radiation for the moment). If a loop is exactly 
circular then it will eventually shrink to a size smaller than r s and form 
a black hole || [7|]. (This is only true for strings that form as a result of 
gauge symmetry breaking. Strings which form as a result of global symmetry 
breaking would radiate all of their energy in the form of Goldstone bososns 



before they shrink to the size of their Schwarzschild radius []19|). 

If the loop is not exactly circular it will still form a black hole provided 
that the size of the perturbations is sufficiently small. Let t s denote the time 
when the unperturbed loop would shrink to the size of its Schwarzschild 
radius. It is clear that if 

\a(r s )y(r s )\ < r s (= a(r s )R s ) (36) 

then a black hole will still form. 

From the analysis of the previous section one could argue that loops of 
string with R << r* could easily form black holes. The argument is that 
since friction exponentially suppresses wiggles on scales smaller than t*, these 
loops would be circular to very good approximation. This would result in 
the copious production of black holes with masses up to 

M ~ fit* ~ A{Gn)~ l m ph (37) 

where m p i is the Planck mass. There are observational upper bounds on the 
density of black holes of masses M~10 14 m p z, which would be evaporating at 
the time of nucleosynthesis or later (see || and references therein). Then, 
from ([T7|), one would be able to put constraints on topologically stable cosmic 
strings of very low tension, 

Gfi < 1(T 14 . (38) 

However, this argument needs to be refined, since for Rq << r* one cannot 
simply use the exponential suppression factor (^) to estimate the effect of 
friction. This is because these loops collapse well before r* and friction has 
not had enough time to damp the perturbations. Therefore, it is important 
to study the behaviour of perturbations on loops with R << t* in more 
detail. 

At very early times, both the loop and the perturbation will be over- 
damped. Neglecting second derivatives in fl3"2| ) we have 

y* L 2 1 



y z R 2 2a(H + h) ' 



(39) 
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Using (p5D we have the interesting relation 



yZ{t)=yZ{to) \^R~) ' (40) 

so the perturbations shrink faster than the radius of the loop, the relative 
perturbation decreasing as the coordinate radius shrinks from its initial value 
Rq. Following similar steps we also have 

y p (*) = y p (<b)(^) i l - ( 41 ) 

To find the limit of applicability of the overdamped approximation, one can 
calculate R and y z from (p5|) and fl39|) and compare them to the damping term 
in equations (|20| ) and (|32|). One readily finds that (|5|) is valid for R « 1, 
while (|40|) is valid for R « L . The time at which R ~ L~ l coincides also 
with the time at which the relevant physical scale comes within the effective 
horizon h~ l and starts oscillating. 

Once a perturbation starts oscillating it is damped very efficiently, so the 
perturbations which will be more difficult to eliminate are those with L = 2, 
which are overdamped up to the time when the loop becomes relativistic, 
R ~ 1. From fl2"5D this happens when R ~ [a/i] -1 , at a time of order Tf, given 
by fl27|). That is, the loop becomes relativistic when 

R 1 /i?o\ 1/3 ; (42) 



R a(T f )h(T f )R 



where in the last step we have restricted attention to loops which are collaps- 
ing in the radiation era. After that, during the relativistic collapse, the loop 
is within the effective horizon and it can be seen that perturbations behave 
very much like they would on a circular loop in flat space. That is, as the 
loop shrinks the amplitude of transverse perturbatins stays constant whereas 
radial perturbations only shrink by a factor of L H. 

As a result, all the supression in the lowest modes L = 2 comes from the 
overdamped regime. From ( |4"T| ) with (R/Rq) given by (f4"2"{), and using that 
r s ~ G[ih~ l (rf) we obtain that at the time t s when the unperturbed loop 
crosses its Schwarzschild radius 

*(«.) ..*<«.>«,„)-. (M 2 ' 3 „ yJM (G , r fM 1/s . (43) 



R s Rq V / -Ro ^ 
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For strings formed at a phase transition, the initial value of the relative 
perturbation y^toj/Ro can be of order one. Then, in order to satisfy (j3G|) 
we need 

However, if this condition is met, the energy of the loop at the moment of 
first collapse is E ~ / ur c ~(G / u) 2 mp;, much smaller than the planck mass, 
which simply means that a black hole will not form. 

Therefore friction by itself is not sufficient to ensure the formation of 
a black hole if we start from an arbitrarily wiggly loop: the lowest modes 
L = 2 are not sufficiently damped. This was not obvious a priory and it is 
in contrast with what would happen if the universe were not expanding. In 



that case h is a constant and from fl4Tj) and (41) all loops whose initial size 
aRo is much larger than (G/i) -1 / 2 /^ 1 would shrink to form black holes of 
mass M ~ 

Even in an expanding universe, some of the loops produced at the phase 
transition might just happen to be circular enough initially that they would 
form black holes, even without friction. With friction included, the number 
of loops that will form black holes is larger. The question is then what 
fraction of the ensemble will go into black holes or, in other words, what is 
the probability of black hole formation Vbh- 

Before we try to answer this question we should know how Vbh is con- 
strained from cosmological observations. Assuming a scale invariant dis- 
tribution of loops with number density at horizon crossing given by [|TJ 
dn(r c ) = vr~ 4 dr c , where v is a parameter of order one or smaller, and given 
that these loops form black holes of mass M = 2n^r c with probability Vbh, 
the number density of black holes is 

MM) = m> (2,G^n h C^-f (if i (£) . (44) 

Here we have included a factor of (r c /t) 3//2 in the distribution to account for 
the dilution of the strings (or black holes) by the expansion of the universe. 
Black holes of mass M ~ 10 10 — 10 n g would evaporate during cosmological 
nucleosynthesis, producing high energy particles which would deplete the 
deuterium and helium when standard nucleosynthesis has almost concluded. 
This process has been studied in Ref . [|T0|| , using a distribution of black holes 
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of the form with the numerical coefficient left as a free parameter. The 
bound that these authors obtained can be translated into 



(27T) 3 ^(G^) 3/2 V bh < 10- 



■2(i 



(45) 



(assuming that the present abundance of deuterium is of cosmo logical origin). 

Note that if Vbh ~ 1> then we would be in conflict with this constraint 
for G/i~10 -18 . It is then important to estimate Vbh- An upper bound can 



be obtained directly from Eq. 
to the probability that 



The probability that y^ts) < R s is equal 



y p 2 (t ) < (G^tJr^Ro. 

Treating 2/2(^0) as a random variable normally distributed, with r.m.s. am- 
plitude of order R , and noting that for L = 2 there are two independent 
modes [the + and the — modes in Eq.fl30l)l, that probability is bounded by 



t x 1 / 3 



-1 2 



4/3 



m. p i 
M 



2/3 



< i(r 9 (G/i) 4/3 , 



(46) 



where in the last step we have used M~10 m p . As a result, the constraint 



flip]) will always be satisfied for cosmic strings in the low energy range 
that we were interested in. 

What can one say about heavier strings? Note that Vbh can actually be 
much lower than the r.h.s. of fl46|) , and the bound fl45|) is likely to be satisfied 
even for large values of the string tension Gfi ~ 10~ 6 . In this case the loops 
that form black holes in the mass range 10 14 — 10 20 m p ; cross the horizon 
much later than so friction only smoothes out the perturbations of very 
large L. This does not help very much when we try to form black holes. The 
smoothness of the strings in this case has to be attributed to other damping 
mechanisms such as gravitational radiation. The probability of black hole 



formation in the absence of friction has been estimated by Hawking and 
others |T2], g], 

Vb h ^(Gfxf. (47) 

In the case of a loop formed by n straight segments, Hawking found p ~ In. 
This form for Vbh results from the fact that in order to form a black hole 
it is necessary to "fine tune" a set of In angles with acuracy given by G\i. 
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Typically, p will be of the order of the number of random variables that 
parametrize the ensemble of loops, since for a black hole to form all the 
parameters have to be fine tuned with an accuracy given by G\x. Polnarev 



and Zembowicz [|T2| studied Vbh for a family of loops containing excitations 
in the first and third harmonics in a Fourier expansion of the solutions to the 
Nambu equations of motion. They studied a two parameter family and they 
found p ~ 2 (with some uncertainty due to the arbitrariness in the definition 
of a probability distribution in the space of parameters.) However, p is likely 
to be larger (and therefore Vbh smaller) since one need not be restricted to a 
two parameter family. In particular, the general solution including the first 
and third harmonics is a five parameter family (see e.g. [|H|), and one may 
expect p ~ 5. 

Other observational constraints come from black holes of mass M ~ 
10 20 m p ;, which would be evaporating at the present time, producing intense 



bursts of 7-rays (see e.g. [[18]] ) . This constraint has been studied in [[12 



whose authors concluded that if p~2 there is conflict with observations for 
strings heavier than GyU~10~ 7 . On the other hand if p > 4 there is no con- 
flict even if G/j, is as large as 10~ 5 . Including friction does not modify these 
results, since for M ~ 10 20 m p / friction is important only for GyU~10~ 20 [see 
Eq.(|57|)]. But for such low values of the tension there is no observable effect 
even if Vbh — 1 ■ 



7 Strings attached to walls 

In this section we consider perturbations on strings which are attached to 
planar domain walls [|IJ. For simplicity, we shall ignore the expansion of the 
universe. Also we shall restrict ourselves to perturbations which lie on the 
plane of the wall. 

Then we are effectively left with a 2+1 dimensional problem, in which the 
string is the boundary between a region of "false vacuum" (the wall) with 
energy per unit area equal to the surface tension cr, and a region of "true 
vacuum" where there is no wall. In the absence of friction, the equations of 



motion are 14 



- ^nx^ = an^. (48) 

Here □ is the covariant d'Alembertian on the worldsheet of the string and 
n M is the space-like unit vector normal to the worldsheet, with n^rin = 1 and 
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n^daX^ = 0. Our sign convention is that n M points towards the wall. It is easy 
to see that equation (|48p has a solution representing a straight string which 
is constantly accelerating due to the tension of the wall. The wall gradually 
disappears as the string moves forward, its energy going into kinetic energy 
of the string ||14|| . The string undergoes the so-called hyperbolic motion, 



asymptotically approaching the speed of light at late times. 

Eq. (48) has to be modified to include the force of friction F v '. To keep 



things general we take 

F" = F(u ± )n u , (49) 

without specifying the dependence in the transverse velocity of the fluid 
u± = u^n^. (Note that u 1 ^ defined in Section 2 is given by u\_ = u±n u .) 
Then the equations of motion read 

_ [l u x v = [a + F(u ± )]n v . (50) 

We can study perturbations on any solution x M (£ a ) of (|50"D using the 



covariant formalism of Ref . [|TJ] , suitably modified to include friction. Since 
only perturbations which are normal to the string are physically observable, 
one need only consider perturbations of the form 

daf = 0(On M . 

Here <fi represents the proper magnitude of the perturbation, i.e. the normal 
displacement as measured by an observer that is moving with the string. 
Multiplying (^) by n v we can write 

-^ a a = (a + F), (51) 

where K a b = n M V a dbX^ is the extrinsic curvature of the worldsheet. Latin 
indices are raised and lowered using the worldsheet metric 

9ab = OaX^dbX^. 

The linearized equation of motion for the perturbations can be found from 
variation of Eq. (pl|) 

pi6K a a = -5F, 

where 5 denotes the variation induced, to linear order in 0, by the small 
perturbation Sx^. Following [14]], we have 5K® = □</> + K ah K a b(\). On the 
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other hand SF = F'(u±)5u±, where F' = dF/du±. Also, 5u_\ 
the change in n M induced by the perturbation is ]14| <5n M = 
Therefore, the resulting equation of motion for is 



- D0 + M 2 <; 
where the "mass" is given by 



u^Sn 11 , and 
■g ab 4>, b d a xv. 



(52) 



M 



K nh K 



ab 



K 



a + F" 



(53) 



In the last step we have used the Gauss-Codazzi relation K ab K ab = (K^) 2 —TZ, 
where 1Z is the intrinsic curvature scalar on the worldsheet. The r.h.s. of eq. 
( |52|) acts like a friction term for the perturbation <ft provided that F' > 0. 
But this condition is always met, since it just means that the force of friction 
increases with the transverse velocity of the fluid. 

In the spirit of Refs. |14|, |15[ , Eq. (|52| ) can be seen as the equation of motion 
for a scalar field <fi which is "living" in the worldsheet of the defect. Now 
this equation has a different mass than in the case without friction, and it 
also has a right hand side in which the field has a derivative coupling to 
some external sources (which are essentially the 'temporal' components of 
the tangent vectors to the string.) 

In the case when the unperturbed string is straight, lying along the y axis 
and moving with trajectory x = x(t), the metric on the worldsheet is 

ds 2 = g ab di a di h = -(1 - x 2 )dt 2 + dy 2 = -dr 2 + dy 2 . 

Here t is the time in the rest frame of the fluid and r is the proper time of 
an observer who is moving with the string. In this case the metric on the 
worldsheet is flat, the curvature scalar vanishes and the effective mass M 2 is 
tachyonic. This is also true in the absence of friction and it essentially 
means that modes with wavelength larger than |M| _1 are unstable. Now, 
including friction, the difference is that the tachyonic mass "switches off" as 
the string approaches its terminal velocity v. Indeed, the terminal velocity 
of the string is determined by the vanishing of the driving force in (|50|), that 
is 

a + F{u ± ) = 0, (54) 
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but this implies that the mass term for the perturbations vanishes. 

When the straight string has reached the terminal velocity v, equation 
( |52"D reduces to 

< 55 » 

where k is the wave number of the perturbation, and j v — (1 — v 2 )~ 1 / 2 is the 
relativistic factor corresponding to the terminal velocity. The friction term 
in the r.h.s. of ([55]) now causes the perturbations to decay exponentially with 
a lifetime which is easily calculable. Taking a force of the form (Q) we have 
F' = (3T 3 . Noting that u± = n^u^ = —vy v , Eq. fl54|) gives 



(3T 3 



In cosmological situations the temperature is always lower than the energy 
scale of the wall (typically of order c 1//3 ), therefore the terminal velocity will 
be relativistic 

v ~ 1 

and the gamma factor will be of order j v ~ aT 3 (ignoring the numerical 
factor 0). Then the r.h.s. of (53) reads 



F' lv a 
fi fj, 

It is then straightforward to see that the perturbations decay with proper 
lifetime given by r ~ ji/a for k >> a / ji and r ~ oir x k~ 2 for k << a j\x. 

A more interesting case is that of a string which is at the boundary of 
a circular hole that has spontaneously nucleated on a metastable domain 



wall (see e.g. Ref.[|l^]). This process is the 2+1 dimensional analogue of the 
formation of true vacuum bubbles in the problem of false vacuum decay. 

Without friction, perturbations on these bubbles (or holes) have been 
studied in ref. [[L4j]. In that case, the unperturbed solution was a circular 



hole whose radius R expands with constant acceleration, 

R 2 -t 2 = (2fi/a) 2 . (56) 

The effective mass for the perturbations was M 2 = —a/2fi 2 , which is tachy- 
onic. Because of the expansion of the hole, any perturbation eventually 
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reaches a wavelength larger than M 1 , at which point it becomes unstable 



and starts growing like oc R ~ t [14|. From an intrinsic point of view, one 



can say that the string is unstable, in the sense that 0, (the perturbation 
measured by a co-moving observer) grows in time. However, an external ob- 
server measures a perturbation which is Lorentz contracted A = 7 _1 0, 
where here 7 = (1 — R 2 )~ 1 ^ 2 . Since from ( |56|) 7 ~ i?, we have A ~ const, at 
late times. The relative perturbation A/i? decreases and the string becomes 
more circular as the hole expands [[U|]. 

Including friction, the main difference will be that the string reaches a 
terminal velocity v, and the Lorentz contraction factor will go to a constant 
at late times. Also, the mass term switches off as the string approaches the 
terminal velocity. Let us see, then, what is the fate of the perturbations 
when friction is included in the dynamics. The string worldsheet is given by 

= (t, Rcos6, Rsin6, 0), and the metric induced on the worldsheet is 

ds 2 = -(1 - R 2 )dt 2 + R 2 d6 2 = -dr 2 + R 2 d6 2 . 

The normal vector is = (1 — R 2 )~ l l 2 (— R, cos#, sin 9, 0). The extrinsic 
curvature can be easily calculated and one finds 

— M 2 = K ab K ab = 7 2 iT 2 + 7 6 (i?) 2 . (57) 

The velocity is bounded by the terminal velocity R < v , so at late times 
R ~ v t and R decays faster than R~ l . Then, neglecting the R term in (|57|) 
the equation of motion for the perturbations takes the form 

d 2 ^ F' lv d<f> L £ 2 ~7^ n ^ 

with (f> = 4>l{t) exp(iL9). 

Note that as the hole expands, the tachyonic mass squared — , j 2 R~ 2 "red- 
shifts" at the same rate as the wave number term L 2 R~ 2 . One might expect 
instability for modes with L < 7^, since for these the total effective mass 
squared is negative. However, at sufficiently late times we can use the over- 
damped approximation and neglect in fl58|). Then a simple integration 
shows that, since R~ 2 decreases faster than r _1 , the amplitude of the pertur- 
bations asymptotically approaches a constant — > const. From an intrinsic 
point of view the asymptotic behaviour of the perturbations is very different 
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from the frictionless case, in which was growing in time. However, now 7„ 
is constant and we have A = 7~V — > const., so from the point of view of 
an external observer which is at rest the behaviour is similar to that of the 
frictionless case. The relative perturbation A/R also goes to zero, and the 
loop becomes more and more circular as the hole expands. 



In this paper we have studied the evolution of cosmic strings taking into 
account the force of friction. The results are conveniently expressed in terms 



For small perturbations around a straight string, a WKB analysis shows 
that they are exponentially suppressed if their wavelength crosses the horizon 
before the time £*. Relative to the frictionless case, the amplitude of the 
perturbations has to be multiplied by the suppression factor (|35|) , which in 
the case of a radiation dominated universe (a = 1) reduces to 



Here A is the physical wavelength of the perturbation at the time and 
7 « 1.2. The suppression of the amplitude of these perturbations is due to the 
oscillation of the perturbations before the time £*. After t* the perturbations 
oscillate with constant physical amplitude. Perturbations whose wavelength 
crosses the horizon after £* are practically unaffected by friction. 

Similarly, for circular loops we should distinguish between large loops with 
r c >> t*, and small loops with r c << t*. Here r c is the radius of the loop 
at horizon crossing. Large loops are unaffected by friction, after they cross 
the horizon they start oscillating with constant physical amplitude r c . Small 
loops do not start collapsing relativistically right after they cross the horizon. 
During the non-relativistic evolution, the velocity at which the string moves 
with respect to the fluid is given by 



where r is the physical radius of the loop, T is the temperature and (3 is the 
numerical factor in (0). Due to the frictional force, these loops are dragged 



8 Conclusions 





(59) 
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by the expansion more than they would in the absence of friction, growing to 
a size larger than r c . However as they collapse they lose energy to friction. 
It is interesting to note that in a radiation dominated universe both effects 
approximately compensate each other, in the sense that at the time when 
they first collapse to a point their energy is given by Ef ~ 2irfir c , just as 
in the frictionless case (in matter dominated universe Ef is actually larger 
than 27r/ir c ). In subsequent oscillations, the loops will keep losing energy, 
up to the time t*. After that they oscillate with constant energy. This 
energy can be obtained from (|28|) in the limit of large times, and one has 
E(t » Q « E / exp(-2^/3r c ) 1 / 3 . 

We have also studied perturbations on circular loops. If the loops are 
larger than i*, then the perturbations behave in much the same way as the 
perturbations on a straight infinite string. Perturbations whose wavelength 
crosses the horizon before i* are suppressed according to (|512D, whereas if it 
crosses later than they are unaffected by friction. If the loops are much 
smaller than £* then it is not correct to use the suppression factor fl59|) , 
essentially because the loops start their relativistic collapse much before t*. 
As a result, the perturbations in the lowest modes are only suppressed as a 
power law in (r c /t*) [see e.g. (f43[)|. 

Because the suppression in the relative amplitude of the lowest modes is 
only power law, the probability that loops of strings become circular enough 
to form black holes is very small (to form a black hole one needs that the 
relative amplitude of the perturbations be of order G/j,, which is very small). 
As a result, the number of black holes of masses M~10 14 mp; (which would 
be evaporating at the time of nucleosynthesis or later) formed by strings 
which cross the horizon before t* is too small to have any observable conse- 
quences (for all values of //). For loops which cross the horizon later than 
t*, friction only eliminates the perturbations of very large wavenumber, so 
the probability that they form black holes is not substantially enhanced by 
friction. 

We have also studied perturbations on strings attached to walls. For this 



we have used the covariant formalism developed in ref. [13], generalizing it 
to include the force of friction. In particular we have studied the case of a 
circular hole which spontaneously nucleates on a metastable domain wall. 
We find that perturbations on the string that is at the boundary of the hole 
are initially unstable, growing at the same rate as the radius of the expanding 
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hole. However, as the loop reaches its terminal velocity v < 1, the instability 
switches off and the perturbations freeze out at a constant amplitude. As 
a result, the relative perturbation decreases in time, and the hole becomes 
increasingly circular as it expands. This behaviour is similar to the behaviour 



that one obtains in the absence of friction |T4j. However the mechanism by 
which perturbations freeze out is very different. In the former case it is due to 
the force of friction, which balances the instability due to the tension of the 
wall. In the latter case the freezing occurred for purely kinematical reasons, 
since the unperturbed string asymptotically approached the speed of light. 
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Appendix A 



Let us solve Eq.([13|) in the WKB approximation. The solution 'under the 
barrier' is given by 



^ = ^2 ex P =t / dr'pi (r' 



T « T k ) (Al) 



a 

± Vpi 

where p\ = [V(r) — k 2 ] 1 / 2 and is the 'classical' turning point 

V(n) = k 2 . 

Since at early times V{r) ~ a 2 r^r 2 (r 4 ,/r) 4 ° ! , it is clear that 
1 = f T k dr' Vl {r') -> _^-^) 2a " 1 -T(k,T k ) + 0(r^), (A2) 
where the function T(k,Tk), which will be specified below, is independent of 

T. 
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Since perturbations are frozen in at early times, the coefficients C± are 
determined by imposing the boundary condition that the co-moving pertur- 
bation y should approach some constant value y as r — > 0. This immediately 
requires 

c + = o, 

C- = y ^/aT^e- r ^\ 

Using the standard connection formulas of the WKB formalism, the so- 
lution in the 'classically allowed' region is 

* = 2y {^P\ ' e~ T ^ cos (£ dr'p 2 (r') - , (r » r fc ), (A3) 

where p 2 = [k 2 — V} 1 / 2 . 

The simplest case, which is also the most interesting, is the radiation 
dominated universe a = 1. In this case the potential assumes a particularly 
simple form V = t 2 t^ 4 , and the turning point is given by r k = (r^k^ 1 ) 1 ^ 2 . 
Then, the exponent T(k, T k ) can be calculated analytically 

r(k,T k ) = -lim 

Integrating by parts we have 

rl (l x | 3 | 

r = A;r fc lim2/ - — = kr k -B(-, -), (A5) 

where B is Euler's Beta function. 

For a 7^ 1, the potential is more complicated and the turning points 
cannot be given explicitly in terms of r* and k. However, for k » r" 1 , the 
potential under the barrier can be approximated by the second term and we 
have 

r k « [ak-W 20 - 1 ] 1 ^. (A6) 

Neglecting the first term in V(r) and following the same steps as before we 
have 

4a -2 V 4a 2 J 
Substituting in (|AT|), ( |A3| ) and using (^) we have the solution (|I~5|). 
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Appendix B 



Here we express the equations of motion for a string with friction in cylin- 
drical coordinates, and derive the equations for perturbations on a circular 
loop. 

The FRW metric is written as 

ds 2 = a 2 (r)(-dT 2 + dp 2 + p 2 d6 2 + dz 2 ). 
The nonvanishing Christoffel symbols are 

a 



r° 

1 00 



90 



zO 



-TP 

1 pO 



a 
a 



r p — —n- r e — n- 1 - F° — -n • 

(X 

The metric on the string worldsheet is given by 

g ab di a de = a 2 (r)[-(l - v 2 )dr 2 + p>d£% 

where 

v 2 = (x) 2 = gijX l x J , 



V 



[x ) Qij% % ^ • 



An overdot denotes derivative with respect to r, and a prime is derivative 
with respect to £. Here §ij is the flat three-dimensional metric in cylindrical 
coordinates, and the indices i,j run over p, 9 and z. 
Substituting in ([I]), the temporal component yields 

e + A{r)v 2 e = 0, 

where A(r) = 2a(H + h) and e = p(l — v 2 )^ 1 ^ 2 . The spatial components 
yield 

1 s j\' r Q/2i 

(Bl) 
(B2) 



p + (l 


-v 2 )A(r)p- 1 - 




y- 


# 2 - 




= o, 


(9 + (1 - 


-v 2 )A(r)6-^- ( 


:9 


' 2 

p 




07/ 
I 5 " 


= 
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= 0. (B3) 

For circular loops we can take z = and the parameter £ proportional 
to the angular variable, £ = Ro6, where Rq is the initial coordinate radius of 
the loop. Then p = R(t), p' = 0, 9 = and we obtain and fl2~T|). 

Let us derive the equations for the perturbations on circular loops. Taking 
z = y z {r6) « R it is immediate from flB3| ) to obtain the equation for small 
transverse perturbations (|32|). The equation for the radial perturbations 
requires more work. First we write the perturbed solution as p = R{t) + A, 
and 9 = (£/Rq) + 6, where A and 5 are small deviations from the unperturbed 
values. Substituting in (pl|) one finds, to linear order in the perturbations 



R + A + (1 - R 2 - 2RA)AR [ 1 + — ) - — + (R + A) — = 0. 
Also to linear order 

_ T 9A 9 R A 1 

[1 + 26') 



2A 2i?A 
1 + ^ + 



R 1-R 2 



where eo is the unperturbed value. Substituting in the previous expression 
and using (|20|) we find 



A - 2^A + (1 - 3R 2 )A(t)A + ^r^A = 0, 
R e z 

which is the equation of motion for radial perturbations (|3l|). 



Figure captions 

• Fig. la Evolution of a circular loop without friction in radiation dom- 
inated universe, as a function of cosmological time t. Here t c is the 
time at which the loop crosses the horizon, r is the physical radius of 
the loop and r c is the radius at horizon crossing. The energy of the 
loop E rapidly approaches the value E c = 2ix^,r c after the loop crosses 
the horizon, and remains approximately constant. We also plot R, the 
velocity of the string with respect to the cosmological fluid. 
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• Fig. lb Evolution of a circular loop with friction in radiation dominated 
universe, with r c = I0~ 3 t*. (Same conventions as in Fig. la.) 

• Fig. 2 Numerical results for the evolution of perturbations on a circular 
loop in radiation dominated universe, ignoring the force of friction. The 
results are plotted as a function of wave number L. The circles denote 
the ratio y p phys ^ L divided by the r.h.s. of equation ( |33"D at the time when 
the loop first collapses. Similarly the crosses denote the ratio of yl hysL 
to the r.h.s. of ([34]). Ignoring the oscillatory behaviour, this ratios are 
of order one, which means that (^) and (0) give a very good estimate. 

• Fig. 3a Numerical results for the evolution of perturbations on a circu- 
lar loop in radiation dominated universe, including the force of friction, 
for R = r*. Circles and crosses denote the same quantities as in Fig. 
3. It is seen that the suppression factor (|3^), depicted as a solid line, 
gives the right answer to very good approximation. 

• Fig. 3b Same quantities as in Fig. 3a, for the case Rq = St*. 
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